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Abstract 

The generation of zonal magnetic fields in laboratory fusion plasmas is predicted by theoretical 
and numerical models and was recently observed experimentally. It is shown that the modification 
of the current density gradient associated with such corrugations can significantly affect the stability 
of the tearing mode. A simple scaling law is derived that predicts the impact of small stationary 
current corrugations on the stability parameter A'. The described destabilization mechanism can 
provide an explanation for the trigger of the Neoclassical Tearing Mode (NTM) in plasmas without 
significant MHD activity. 

PACS numbers: 
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I. INTRODUCTION 



In nature as well as in experimental devices, plasma turbulence generates coherent struc- 
tures, which can significantly affect the overall behavior of the system. An example in 
magnetic fusion plasmas is the occurrence of zonal electric fields or zonal flows (the velocity 
is generated through E x B drifts). The zonal flows are predicted by electrostatic turbulence 
theory [l] and are observed as mesoscale oscillations in the velocity field, which fluctuates 
in the radial direction while being homogeneous in the toroidal and poloidal direction (i.e. 
m = 0, n = 0, where m and n are the toroidal and poloidal wave numbers in a toroidal 
geometry). Such structures are peculiar to turbulent systems since they are linearly stable 
and can only exist as a result of nonlinear wave interaction. 

The importance of the zonal flows as a self-regulating mechanism for the turbulence 
and^plasma transport is widely accepted. Only recently 
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y, however, experimental observations 
, 3,1s] have highlighted that turbulence 



and nonlinear electromagnetic simulations Ref. [ 
in fusion plasmas can also generate, together with the well documented zonal flows, zonal 
fields. These, similar to the zonal flows, are axisymmetric sheared band-like structures in 
the magnetic field. The zonal fields were predicted by several theoretical works, a detailed 
review of which can be found in Ref.jlj]. In these works, zonal magnetic fields are shown 
to arise from finite (3 (the ratio between the thermal and the magnetic pressure) drift-wave 
turbulence when electron inertia effects are included. 

The radial oscillation of the zonal magnetic field is associated with zonal currents flowing 
in the direction of the confining magnetic field. These turbulence generated current cor- 
rugations perturb the m = 0, n = component of the current density and its gradient, 
therefore affecting the stability of the modes driven by it, such as the tearing mode [7|. The 
reconnection of the magnetic flux caused by the tearing instability leads to the formation of 
so-called magnetic islands which, if macroscopic, can significantly affect the performance of 
fusion experimental device, as their particular topology enhances the radial transport and 
therefore reduces the confinement. 

The stability of the tearing mode is conveniently measured by the stability parameter 
A', a positive value of which corresponds, in the simplest theoretical models, to an unstable 
mode [7, s|. When pressure effects are neglected and the boundary conditions of the problem 
are fixed, A' is a function of the m = 0, n = component of the current density and of the 
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wave vector o 
In Refs. 



the perturbation, k, only. 
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ll| it was shown that, even if A' depends on the global features of the 



current density profile, it is strongly affected by the local structure of the current density at 
the resonant surface where k-B = 0. This observation suggested the possibility of controlling 
the tearing mode by applying small localized current perturbations around the reconnecting 
surface. In a similar way, however, turbulence generated current corrugations can provide a 



source of free energy for the tearing mode and modify its stability [12] . As a consequence, 
the zonal fields could be a possible triggering mechanism for Neoclassical Tearing Modes, as 
they could push the island above the seed limit. Implied in this reasoning is the assumption 
that the corrugations are sufficiently coherent for the tearing mode to have time to respond. 

In this work we describe the effect of current corrugations of given shape and amplitude 
on the calculation of A'. By using a simple slab model for the zonal fields we obtain 
an analytical scaling for the variation of the stability parameter, 5 A', with respect to the 
amplitude and the wave length of the "zonal" flows that generate the current corrugations. 
We then show numerically that the scaling applies also to more relevant cylindrical cases. 
We find that even relatively small current corrugations can induce large modifications of A' 
and therefore strongly affect the stability of the tearing mode. 

A limitation of our work is the assumption, done for sake of simplicity, that the pressure 
is flat in the island region. Pressure gradients at the reconnecting surface, however, can 
play an important role in the definition of the stability parameter, as shown in Refs. [3, 
We will address this subject, together with a discussion of the generation of the "zonal" 
pressure, in a future publication. 

The paper is organized as follows. In Section II we describe the basic model behind 
the calculation of A' for a given equilibrium current density. We then describe the effect of 
small modifications of the current density on the stability parameter by using a perturbative 
approach. These results are applied to a slab configuration in Section III, where a simple 
analytic scaling for 5 A' is found. In Section IV we apply the scaling to a cylindrical problem 
and we verify its validity using a numerical code. Our results are used to investigate the 
onset of the NTMs in the presence of current corrugations in Section V. Finally, in Section 
VI we draw our conclusions. 
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II. MODEL 



The 2D fluid equations employed in our analysis are well suited to describe plasmas 
confined by a large toroidal magnetic field. We restrict our attention to configurations with 
small inverse aspect ratio, e = a/Ro, where a and Ro are the minor and major radius of 
the machine. This leads to a simplification of the geometry and, in particular, it allows the 
use of a cylindrical tokamak approximation. In the following, the cylindrical coordinates 
(r,6,z) represent the radial, poloidal and "toroidal" direction, respectively. Furthermore, 
we order f3 as e, which implies that the electromagnetic effects are relatively small, although 
not negligible. 

We represent the magnetic field as: 

B = B z e z + Vtpp x e 2 . (1) 

where B z = const is the component of the magnetic field in the "toroidal" direction (i.e. 
along the cylinder axis) and ip p is the poloidal magnetic flux. In order to describe a single 
helicity tearing mode it is convenient to introduce an helical coordinate: £ = 9 — ^j^z. 
In the same way, we express the magnetic flux in terms of its helicoidal component, which 
is given by: ip — ip p + ^jfo ■ Finally, the helicoidal magnetic flux can be decomposed in 
its equilibrium and perturbed part (the tearing mode) and the latter can be expressed by a 
Fourier series: ip = ip e q(r) + '^^ip m (r)e l ^' yt ~ m ^ . The magnetic flux and the "toroidal" current 

m 

density, J, are related through Ampere's law, so that (An/c)J = — Vj_Vp = — + 7^"' 
where Vj_^ = r~ 1 d r (rd r ip) + r~ 2 d 2 ip and c is the speed of light. The mode resonant surface, 
where the flux reconnects, is such that the magnetic winding number (the safety factor), 
q{r) = is equal to m/n (Be is the poloidal magnetic field). The safety factor then 

can be easily expressed as a function of the equilibrium "toroidal" current, (An/c)J eq = 
-r^K^ + H: 

= c r 2 B z 
q 47rR jdr[rJ eq (r)Y 1 ) 

In a high temperature plasma with low resistivity, and within our approximations, the 

structure of the linear tearing mode far from the resonant surface can be obtained by solving 

the scalar equation e z ■ V x (J x B) = (i.e. the projection of the curl of the plasma 

momentum balance without inertial or viscous effects). Using Eq{T] [14] and the definition 
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of the helicoidal flux it is possible to cast this equation in the following form: 

[J,V] = 0, (3) 
where the operator [A, B] = r~ l (d r AdeB — dgAd r B). The linear version of Eqj3]is: 

= 0, (4) 




where the helicoidal magnetic field is given by: — ip' eq = — the prime represents 

derivation with respect to r, the subscript m is dropped and the factor An/c is absorbed in 
the definition of the current. It is clear now that the ideal MHD approximation employed 
so far does not hold in the neighborhood of the resonant surface, as the term proportional 
to 1/ip' eq in EqH] yields a singularity. 

The singular behavior is regularized by allowing for a small but finite resistivity, which 
smooths out the current divergence within a narrow boundary layer centered around the 
resonant surface. As a consequence, we can identify two separated regions, where different 
simplifications of the model equations apply: an "outer" linear ideal region, where EqJD 
holds, and an " inner" linear or nonlinear resistive boundary layer, the narrowness of which 
allows geometrical simplifications, although non-ideal physics must be retained. The solu- 
tions of the equations in the two regions must match in a so called overlapping region. In 
other words, the "outer" solution provides the boundary conditions for the "inner" solution. 

To simplify the matching procedure, it is convenient to introduce the stability parameter, 

A': 

dip 
dr 



A' 




(5) 



r=r s +e 

where ip Ts is the value of the perturbed magnetic flux at the resonant surface. Clearly, A' 
is fully determined once the outer magnetic flux is found by solving EqJU and is a function 
of m, of J' eq /ip' eq and of the boundary conditions imposed on ■0. The stability parameter 
defines the stability of the linear tearing mode since, generally, a positive A' corresponds to 
an unstable mode Q]. Similarl y, a lso in non- linear theory, A' can be considered a measure 



of the stability of the mode 



id- 



We focus now on the modifications to the stability parameter due to local changes of 
the profile of J' eq /ip' eq - We first assume a given equilibrium magnetic flux ipeqo, associated 
with current density J eq0 (and therefore to a safety factor q ). We then superimpose a 



"zonal" field, Sip eq , which can be thought as the m = 0, n = component of the magnetic 
flux perturbation (and is accompanied by a "zonal" current SJ eq and a Sq). The physical 
origin of the corrugations, as discussed in the introduction, is suggested to arise due to a 
combination of short scale turbulence and zonal flows. We again emphasize that, in this 
calculation, we assume that they are sufficiently consistent in time and space such that a 
new "equilibrium" is temporarily created (hence the subscript eg). Obviously, there needs 
to be a constant source of corrugations or otherwise they would dissipate on a time scale 
proportional to a 2 ft], where a (of the order of the collisionless skin depth) is the characteristic 
length scale of the zonal field. In other words, Sip eq directly modifies ip eq o, thus generating 
the new "equilibrium": i[) eq = ip eq0 + Sip eq (together with J eq = J eq0 + 5J eq and q = q + Sq). 
This new configuration affects the perturbed magnetic flux through EqJH and therefore the 
stability parameter through Eqj5l In order to explicitly separate the effect of the "zonal" 
field, we write the solution of EqJUas ip — V'o + <^0j where ip is the eigenfunction associated 
with the old equilibrium, while Sip is the net effect of the corrugations on the mode. Similarly, 
also the stability parameter becomes: A' = An + 5 A'. 

n n 

Following Refs.[9|, UQfl, we assume that the localized current density corrugation is such 
that Sip <C V'O) while Sip' ~ ip' . Furthermore, we assume that Sq <C go, which is valid when 
J dr(rSJ eq ) f dr(rJ eq0 ) (as in the case of a localized zero-average oscillating SJ eq ). With 
these approximations [i.e. neglecting terms O(Sifj/ip , Sq/q )] we obtain from EqJU 

Tr ( r ^) = - r ^ ' (6) 



and consequently, after integration: 

d 



dr 

±$ 

dr 



rs-e SV ~ 

dr-^—^iPo, (7) 
, dT —^^ (8) 



r 



where we have employed the fact that any change in the perturbed flux function must vanish 
far from the region where the current corrugation is located (i.e. dSip/dr = both for r = 
and for r = oo). From EqsEl and [8] we obtain: 



5 A' = — lim 




i P r rdr s Jk$L 

r s Jo V ^9 As 



(9) 
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where P stands for the principal value of the integral (necessary because the integrand is 
singular around the resonant surface). 

In the next Section we explicitly calculate 5 A' for a given current density corrugation in 
slab geometry. This allows us to obtain a semi- analytic scaling of 5 A' with respect to the 
features of 5ip eq . 



III. SLAB CASE AND 6 A' SCALING 



In order to understand how a current corrugation affects the stability parameter, we first 
investigate its effect in a simple slab configuration. The problem is analyzed in a double 
periodic rectangular box, where the normalized "radial" variable is — tt < x < tt, and the 
normalized "helicoidal" variable is —Tr/k y < y < ir/k y , where k y is the aspect ratio of the 
box (and also the wave number of the mode). We assume a given normalized equilibrium 
current density, J eq o = cos(x) [consequently, ipeqo = cos(x)]. With this choice, resonant 
surfaces form in the center of the box at x — 0, and at the edge at x = ±7r, where ip' eq0 = 0. 

In this Section we study an elementary current corrugation, 5J eq = AK 2 cos(Kx + a), 
generated by a "zonal" magnetic field [1] of given shape: 5tjj eq = Acos(Kx + a). Here A 
is the (normalized) amplitude of the "zonal" field, K its wavelength, and a its phase with 
respect to the resonant surface. In order to investigate its effect on the central resonance, 
we reduce Eqj9] to a form that is appropriate to slab geometry: 

SJ 'eg "00 



5A' 



dx 



AeqO VV 



(10) 



Note that x is normalized to a typical length scale and so is 5 A'. Furthermore, slab geometry 
requires symmetric boundaries far from the reconnecting surface, so that the lower extreme 
in Eqj9], r = 0, is replaced by x = — oo in EqJTQl 

For the equilibrium and the corrugation studied here, it is straightforward to obtain: 



5 A' = AK 3 I cos(a) 



tt/2 



dx 



■k/2 



sin(i^x) ipQ 



sin £ 



+ sin(a)P 



tt/2 



dx 



n/2 



cos(Kx) TpQ 



sin (x 



;n) 



As a consequence of the parity of the integrand, the second term on the right hand side of 
the previous equation is equals to zero. At the same time, the first term on the right hand 
side does not contain any singularity at the resonance, which implies that the principal value 
notation can be dropped. Note also that the integration limits are reduced to x = ±7r/2, 



which is halfway between the edge and the central resonance but still asymptotically far 
from the resistive layer (supposed to be infmitesimally small). We restrict our attention 
to tearing modes close to marginal stability, for which the constant- 1 ^ approximation Q is 
valid, i.e. ip' Q <C ipo in the vicinity of the resonant surface. As a consequence, for these modes 
we can assume that ipo(x) /ip rs = 1. Therefore, the integral f^, 2 dx[sm(Kx)/ sin(x)] can be 
calculated for different values of K. For K > 5 its value is roughly constant and can be 
approximated with it (the maximum error in this range is around 10%, and it reduces for 
large K). We remark that the matching theory described here is valid if 1/L -C K < 1/A, 
where L is a macroscopic length scale and A is the resistive boundary layer width. Indeed, if 
KX is too large the scale of the corrugation is of the same order of the boundary layer, thus 
making impossible to employ a perturbative technique, while if it is too small the effect of 
the corrugations is evanescent. 

To summarize, the modification of the stability parameter of a constant-^ tearing mode 
depends linearly on the amplitude of the "zonal" magnetic field and has a strong cubic 
dependence on its wavelength : 

5 A' « Ccos(a)AK 3 , (12) 

where C is a constant depending on the equilibrium and the system configuration (in the 
simple case treated here C — 7r) . When considering the features of the current corrugation, 
the scaling is linear both in its wavelength (K) and amplitude (AK 2 ). Although this scaling 
was obtained with a basic model, it applies also to more complicated geometries, as will be 
shown in Section IV. 

As an application, which is also useful to verify the validity of EqJT2l we study numeri- 
cally the dispersion relation of a tearing mode generated by a corrugated equilibrium. The 
calculation is performed in the same 2D slab double periodic box, equilibrium and corru- 
gation described above. The equation that we solve are the plasma vorticity equation and 
Ohm's law (resistive reduced MHD model): 

dU „ . 

* = V » J - (13) 

^ = s-'M, - J), (14) 

where all the quantities are normalized: the lengths with respect to typical macroscopic 
scale, the time to the relevant Alfven time, S is the Lundquist number, and all the other 

8 



variables accordingly. The operator d ■ ■ ■ jdt = d ■ ■ ■ jdt + Ve x b • V • • • contains the E x B 
drift velocity, Vy = B/|B| ■ V is the parallel gradient, and the vorticity U — e z ■ V x Ve x b- 
The current density is related to the magnetic flux through Ampere's law, which in a slab 
gives: J = —V 2 ip. 



The theoretical dispersion relation appropriate for this model is 191 ]: 



7 s/4 _ 7 V4 65 --i = .48(A' + 8M)S-V A k l J\ (15) 

where b = [J" q o(0) + 5J" g (0)]/[J eg o(0) + 5J eq (0)] (note that for b = we recover the classic 
solution of Ref.p]], called FKR solution in the following). Inspection of EqJT5l shows that 
current corrugation affects the growth rate through both 5 A' and b. 

The linearized version of the system fT3lfT4l is solved using a finite difference eigenvalue 
numerical code, benchmarked with analytical cases. In all the simulations, the Lundquist 
number is S — 10 3 . In absence of corrugations, the stability parameter is uniquely defined 



by k y according to the relation: A' Q = 2(1 - k 2 y ) 1 / 2 tan[(l - fcJ) 1/a 7r/2] (see e.g. Ref.[2(|). We 
start by choosing k y = 0.94, so that the "smooth" equilibrium is characterized by Aq = 0.41 
and it is therefore tearing mode unstable. Solving EqJT5l we find that the perturbation has 
growth rate 7 = 0.0033 (note that the growth rate is normalized with respect to the Alfven 
frequency) and the resistive layer has a width A = 0.045. 

The presence of a corrugation with wavelength K can significantly change the growth 
rate of the mode, as shown in Fig{TJ In the figure the crosses, circles and squares represent 
the numerical solutions obtained with the linear code for corrugations of amplitude 5 x 10~ 5 , 
1 x 10~ 4 , 2 x 10 -4 , respectively. When K < 40 (i.e. KX < 1.82), these values lead to current 
corrugations of maximum amplitude smaller than 7%, 15% and 30% of the equilibrium 
current density. We remark that even small amplitude corrugations can strongly modify the 
growth rate of the mode. For example, in the case shown here a corrugation that amounts to 
4% of the amplitude of the equilibrium current density (e.g. A = 1 x 10~ 4 and KX = 0.91) 
can lead to a growth rate three times as large than in the "smooth" case. For values of 
KX smaller than 1 (K < 20 in our example), the theoretical predictions obtained from the 
dispersion relation Eq{15] together with the scaling EqJT2l are in excellent agreement with 
the numerical data. The disagreement for larger values of the corrugation's wavelength is 
due to the limitations of the validity of the scaling. Furthermore, also the dispersion relation 
Eqjl5]is not correct when the constant-'?/' approximation does not hold anymore. 
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FIG. 1: (Color online). Growth rate of the tearing mode as a function of the corrugation wavelength 
K. The solution is obtained for S = 10 3 , A = 0.41, A = 0.045 and a = 0. Cr osses, circles and 
squares represents cases with A equals to 5 x 10~ 5 , 1 x 10~ 4 , 2 x 10 -4 respectively. The solid lines 
show the theoretical prediction for the different cases. The thin vertical line represents the limit 
of validity of the theoretical prediction (KX = 1). 

To complement this analysis, we have independently calculated with a different numerical 
code the value of 8 A' associated to the corrugations and the equilibrium described above. 
The code, that is also used to obtain the cylindrical results presented in the next Section, 
can solve EqJH or its slab version by using a "shooting" algorithm. As Figj2] shows, the 
numerical values of 8 A' perfectly match the cubic behavior predicted by the scaling in 
EqJT2~| thus confirming its validity. 

We complete this study by showing that a current corrugation with the appropriate phase 
can drive a linear tearing mode even when the "smooth" equilibrium would assure stability. 
In Figj3]we describe three cases, with A' = (red squares), Aq = —0.29 (blue circles) and 
A' Q = —0.8 (green crosses), for a corrugation of amplitude A = 1 x 10~ 4 . As expected, if the 
wavelength of the corrugation exceeds a critical threshold, depending on A' , the instability 
can be excited. We conclude remarking that, for a = 0, the choice of a positive A leads to 
destabilization, while a negative value would produce the opposite effect. 
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FIG. 2: (Color online). 5 A' as a function of the corrugation wavelength K, for A = 0.41 and 
a = 0. Crosses, circles and squares represents cases with A equals to 5 x 10 -5 , 1 x 10 -4 , 2 x 10 -4 
respectively. The solid lines show the theoretical prediction (Eq JT2"|) for the different cases. 

IV. EFFECT IN CYLINDRICAL GEOMETRY 

The analytic results in Section III suggest the possibility that the stability of the tearing 
mode can be significantly affected by the presence of small scale, small amplitude current 
corrugations. The scaling that we have obtained is verified by numerical codes in slab 
geometry, and its validity is confirmed in this Section for cylindrical geometry, which is 
more relevant for experimental applications. In order to proceed, we perform a numerical 
study with the "shooting" code introduced in the previous Section. 



2l| given by: J eq0 = J /[l + 



We assume a "rounded" equilibrium current density 
(r/ro) 4 ] 3//2 , which implies an equilibrium magnetic flux: ip eg o = — 0.25JorQ sinh _1 [(r/ro^ 



21 



Here Jq represents the current density at the cylinder axis and r$ is a measure of the width 



of the current channel. In the following we take Jo = 2/0.7 and r = v/1/5, so that the 
mode m = 2, n = 1 resonates at r s = 0.73 and m = 3, n = 2 at r s = 0.61. In our 
model, the "zonal" field that generates the current corrugation is sinusoidal and localized 
around a radius R by a Gaussian envelope of width p (taken as 0.2 in all the simulations): 
5^ eq = Acos[K(r - R^e^-^/p 2 . 

We investigate the effect of the corrugation on m = 2, n = 1 and m = 3, n = 2 
modes, which are stable in the "smooth" configuration since A' 21 = —2.19 and A' 3 2 = 
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FIG. 3: (Color online). Growth rate of the tearing mode as a function of the corrugation wavelength 
K. The solution is obtained for S = 10 3 , A = 1 X 10~ 4 and a = 0. Crosses, circles and squares 
represents cases with Aq equals to —0.8, —0.29 and respectively. The solid lines show the 
theoretical prediction for the different cases. 

—3.28. We observe that the effect produced by 6ip eq on the stability parameter changes 
as the localization radius of the corrugation is modified and it reaches its maximum when 
R = r s (see FigsJH and [5]). This is expected, since the rigid shift around the resonance of 
the sinusoidal perturbation is equivalent to a change in the phase a in in the slab model. 
Furthermore, the comparison between the solid (A = 5 x 10~ 5 ) and dashed lines (A = 
2.5 x 10~ 5 ) in FigsJH and [5] confirm also in cylindrical geometry the linear dependence of 
5 A' on the perturbation amplitude, as given in EqJT2l Similarly, also the scaling in the cube 
of the corrugation wavelength is verified by measuring the maximum amplitude of the solid 
(K = 40) and the dash-dot (K = 30) lines. 

Finally, it is interesting to evaluate the constant C for the cylindrical equilibrium investi- 
gated. In order to do that, we divide the maximum value of 5 A' (obtained at R = r s ) by the 
amplitude of the corrugation and the cube of its wavelength. With this procedure we find 
that C = 2.5 for m = 2, n = 1 and C =1.7 for m = 3, n = 2. While both values are pretty 
close to the slab estimate, this result suggest a dependence of the effect of the corrugation 
on the wave-vector of the driven perturbation. This dependence is easily explained by the 
presence of the n/m term in the definition of ijj' eq0 that, in cylindrical geometry, appears in 
the integral of EqJTOl 
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FIG. 4: (Color online). A' as a function of the localization radius R, for the mode m = 2, n = 1 
and the equilibrium described in the text (for which A' 2 1 = —2.19). The solid, dashed and dash-dot 
lines describe cases with A = 5 x 1CT 5 and K = 40, A = 2.5 x 10~ 5 and K = 40, and A = 5 x 10~ 5 
and K = 30, respectively. 
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FIG. 5: (Color online). A' as a function of the localization radius R, for the mode m = 3, n = 2 
and the equilibrium described in the text (for which A3 2 = —3.28). The solid, dashed and dash-dot 
lines describe cases with A = 5 x 10" 5 and K = 40, A = 2.5 x 10~ 5 and K = 40, and A = 5 x 10~ 5 
and K = 30, respectively. 
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V. NTM TRIGGERING MECHANISM 



In this Section, we analyze how the zonal fields affect the onset of the Neoclassical Tearing 
Modes. In order to do that, we employ a reduced version of the generalized Rutherford 



equation 



22 



23| : 



dw . . . . . _ / w 
t v —kA' + 5A' + C — 
at \ w 2 + w 



WW 



pol 



/,„./ + W pol 



(16) 



Consistently with Section III, the island width and the stability parameter are normalized 
to a macroscopic length scale (e.g. the resonant radius). For sake of simplicity, we retain 
in the model only the fundamental terms responsible for the stabilization/destabilization of 
the NTM. In particular, we neglect the current shape term, tha t strongly affects the island 
saturation but has little effect on the seed island problem [24], 



curvature effects 



25] and the stabilizing toroidal 



26] . Furthermore, we do not take into account the rotation of the magnetic 



island in a self-consistent manner 



27 



281 ]. Indeed, in EqJT6]neither the bootstrap term (third 



on the RHS), nor the polarization term (fourth on the RHS) are explicit functions of the 
island rotation. The coefficient w^w, is a measure of the reduction of the drive associated 
with the perpendicular transport [29] , w po i is the cut-off due to the banana orbit effect [3^] 
and Co is a constant proportional to the poloidal (3 [22]. 

We remark that our purpose is to sketch with an heuristic model the mechanism that 
could allow NTM formation. A detailed analysis would require a self consistent treatment 
of Rutherford equation coupled with an equation for the rotation frequency of the magnetic 
island. The theoretical uncertainties on the exact form of this coupling make an attempt to 
employ more complicated models to study the problem at hand a futile exercise. 

In absence of zonal fields and if Aq < 0, Eq{16] predicts a minimum island width, w see d, 
below which the perturbation is always stabilized and the NTM does not reach macroscopic 
size. Conversely, when the threshold is exceeded, the dynamical evolution takes the island to 
its stationary saturation amplitude, w sa t- To destabilize the NTM an island of sufficient size 
is therefore required. This "seed" island could be produced by some other MHD activity in 
the plasma, such as sawteeth, fishbone instabilities or Edge Localized Modes. In the following 
we investigate NTM dynamics without a seed [i.e. w(0) = 0] island but in presence of zonal 
fields. 

From experimental observations 

aa 

, we expect that the zonal fields maintain their 
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FIG. 6: (Color online). NTM triggering mechanism, (a) Approximate boundaries of the relevant 
solutions of EqJ16t obtained from the numerical solution of Eqll6l Parameters 5 A' and £ in zone 
1 do not produce any island, in zone 2 yield to small (w < w see d) oscillating islands, while in zone 
3 islands reach full NTM saturation. Physical parameters used described in the text, (b) Time 
evolution and (c) phase diagram of an island with parameters in zone 2 (£ = 0.05, 5A' max = 2.94, 
solid line) and in zone 3 = 0.1, 5A' max = 2.94, dashed line). The thin solid lines in (c) show 
Trjdw/dt for the maximum and minimum value of 5A' , while the dash-dot line for 5A' = 0. 

coherence for a limited amount of time, r z f. This fact is modelled by assuming a temporal 
variation of the modified stability parameter, so that 5 A' = SA' max cos[(2tt /r z f)t], where 
5A' max is a constant. A critical parameter of the problem is therefore the ratio between this 
time scale and that typical of the evolution of the NTM. From EqJT6]we infer that this time is 
a fraction of the local resistive time (calculated with the neoclassical resistivity and the minor 
radius) t^tm ~ O.lr^, since the saturated island width can be usually estimated around one 
tenth of r s . We expect that rapidly varying zonal fields, such that £ = r z f/r n -C 1, will not 
be an efficient trigger for the NTM as their effect is averaged to zero. On the other hand, 
if T z f/r v is of the order or greater than one (although the latter possibility is difficult to 
achieve in fusion plasmas) the zonal fields can drive w above its seed island limit. 

The solution of EqJT6] leads to three basic dynamic behaviors for a stable configuration 
(Aq < 0) in the presence of a corrugation. When SA' max < \A'\ the corrugation is not strong 
enough to produce any island (region 1 in FigJH](a)). If the modified stability parameter is 
sufficiently large, but smaller than a critical value <5A^"^, or r z f/r v < ( crit , an oscillating 
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island appears (region 2 in FigJB^a)). Its maximum width remains below w see d and is roughly 
proportional to r z /. Finally, when the critical values are exceeded, the full NTM is triggered 
and the system settles in a state with w oscillating around w sat (region 3 in Figlf^a)). The 
approximate boundaries of these regimes [obtained solving Eq{TB] numerically] are sketched 



in Figl^a), for a case with A' : 
reasonable experimental values 
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., Cp = 0.25, Wb 00 t = Wp i = 0.025, which correspond to 



321 ] . For clarity, in Fig 16] (b) we show the time evolution 



of two cases with 5A' max = 2.94, one just below ( crit ( ( = 0.05 for the solid line) and one 
above (( = 0.1 for the dashed line). For the same two solutions, we plot in Figj6] (c) the 
phase diagram, which describes how dw/dt and w are related during the evolution of the 
perturbation. Two dots in FigJ^a) display the position of these cases in the SA' max — ( 
space. 

Obviously, in a realistic situation the time evolution of 5 A' would be erratic and Eqfl6l 
would become a stochastic differential equation. We expect that this would make the trigger- 
ing mechanism easier, with a lower threshold for both the critical ( and 5 A'. An interesting 
discussion of the stochastic behavior of Rutherford equation was given by Itoh et al. in 
Ref. 23j. Furthermore, the calculation presented here is conservative and the critical values 
obtained are very stringent, since we have assumed an initial condition for Eq{16]with no 
magnetic island. In reality, the presence of a small, but finite seed island island [w(0) < w see d] 
would move the stability boundaries in Fig® broadening the unstable region 3. 



VI. CONCLUSIONS 



We have investigated the effect of current corrugations on the stability of the tearing 
mode. We have obtained a theoretical scaling which relates the change of A' to the charac- 
teristics of the corrugation, such as its amplitude, typical scale length and phase with respect 
to the resonant surface. The scaling has been tested on a simple slab linear problem, giving 
excellent agreement between theoretical predictions and numerical data. Then, by using a 
shooting code, we have investigated the effect of the corrugations in cylindrical geometry, 
and found qualitative and quantitative agreement with the analytic scaling. Finally, we have 
addressed the problem of the seedless trigger of a NTM through a current corrugation. 

From our results we conclude that even small modifications of the original equilibrium can 
lead to significantly different values of the stability parameter and can therefore strongly 
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affect the onset and the dynamics of the mode. For example, a current corrugation of 
amplitude around 5% of the equilibrium and with scale length comparable with the resistive 
layer can triple the growth rate of the tearing instability. The study presented in this paper 
has two clear implications for tokamak experiments. 

First, it implies that the calculation of the stability parameters from experimental data 
is an extremely delicate procedure. Indeed, the current density profiles used to evaluate A' 
cannot be directly measured. While the magnetic field can be deduced from Motional Stark 
Effect measurements, the best spatial resolution achievable is around 5cm, too coarse to 
rule out the possibility of current corrugations. Another option is to reconstruct the current 
profile from temperature and density data by using equilibrium codes. However, the error 
bars of the measurement and the approximations used in the current reconstruction make 
it virtually impossible to obtain a precise value of the stability parameter. 

The second implication is related to the trigger of the NTMs in a relatively quiescent 



plasma. Experiments in several machines 
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32 



33] have reported Neoclassical Tearing 



Modes growing without strong MHD precursors (sawteeth, fishbones or ELMs) which could 
generate seed islands. The theory presented in this paper suggests the possibility of NTMs 
triggered by the microscopic turbulence via the generation of slowly evolving (on the NTM 
time scale) current corrugations. The corrugated configuration could be linearly unstable 
and generate a seed island, which would then be sustained by Neoclassical physics. At this 
point, the pressure flattening in the island region could remove the drive for the turbulence 
and therefore for the zonal field, which would decay, leaving the island at the NTM saturated 
level. The heuristic model that we have employed shows that a seedless NTM can indeed 
be excited if the zonal flow that produces the corrugation is sufficiently strong and slowly 
evolving. Furthermore, even when the zonal field is not able to directly drive the NTM, it 
can significantly reduce the seed island width (i.e. the plasma becomes more sensitive to 
the presence of small islands). 

It is useful here to give some experimental values for the parameters of our NTM trigger 
theory. For a typical tokamak plasma, the resistive time is of the order of a few second, 
hence the islands saturate in roughly 50-100 ms, while the peak frequency for the plasma 
turbulence is around 100 KHz (as observed for example in Ref. 34|). According to experi- 



mental observation of zonal fields by Fujizawa et al. [3j, current corrugations evolve roughly 
50 times slower than the drift-wave turbulence. This is in agreement with the theoretical 
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prediction that zonal fields are mesoscale objects. The extrapolation of these results to typ- 
ical tokamak plasmas leads to an expected coherence time for the corrugations around 1-2 
ms. As a consequence, in standard conditions, the value of ( is in general below the critical 
value predicted with our model, which is in agreement with the fact that, in tokamaks, seed- 
less NTMs are the exception rather than the rule. However, current corrugation of sufficient 
amplitude occurring in particular plasmas (e.g. where the local resistivity is particularly 



high or the zonal fields are slow) could explain the results in Refs. 3l|, |32|, |33|. Furthermore, 
assuming that the zonal field time scale will be roughly unchanged in the next generation 
tokamaks, such as ITER, our results seem to suggest that these machines should be stable 
to spontaneous NTMs, as their resistive time is of the order of hundreds of seconds. 

A limitation of this work is the assumption that the pressure profile (and therefore the 
pressure perturbations) does not play any role in the calculation of A'. A "zonal" pressure 
is generated by drift-waves when Finite Larmor Radius effects are considered and should 
become significant when ITG instabilities are present. Pressure gradients might be relevant 
as they contribute to Eq|4] with a term that is inversely proportional to {ip' eq ) 2 , and are 
therefore very important in the neighborhood of the resonant surface. As pointed out in 
Refs. 3_5| this could make even harder to properly evaluate the stability parameter. The 
inclusion of this effect in our calculation will be discussed in a future publication. 
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